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ABSTRACT 

form  solution.  In  most  of  the 
quantum  mechanics  etc.,  the 
not  available,  we  go  in  for 

numerical  solution  of  the  partial  differential  equations  by  various  methods. 
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1.1  INTRODUCTION 


Many  partial  differential  equations  cannot  be  solved  by  analytical  methods  in  closed 
research  work  in  fields  like,  applied  elasticity,  theory  of  plates  and  shells,  hydrodynamics, 
research  problems  reduce  to  partial  differential  equations.  Since  analytical  solutions  are 


The  general  linear  partial  differential  equation  of  the  second  order  in  two  independent  variables  is  of  the  form: 

A (x'y&+  B c@r)0  + F = 0 


Such  partial  equation  is  said  to  be: 

(l)Elliptic  if  — 4 AC  < 0,  (2)  Parabolic  if  E 2 — 4 jIC  = 0, 

(3)Hyperbolic  if  — 4 AC  > 0, 

The  one  dimensional  heat  equation,  namely, 

— = a — — Where,  a — — is  an  example  of  parabolic  equation 
St  pc 

2 1 

Ct  — - , the  equation  ( 1 .20),  becomes,  — — a — = Q, 
c &x 1 d r 

Herej4  = 1,  B = 0,  C = 0, 

B~  — AAC  = 0.  Therefore,  it  is  parabolic  in  all  points. 

The  solution  of  this  equation  is  a temperature  function  It  (:t,  £j  which  is  defined  for  values  of  X from  O to  / and 
for  values  of  time  t from  0 to  tfi.  The  solution  is  not  defined  in  a closed  domain  but  advances  in  an  open-ended  region  from 
initial  values,  satisfying  the  prescribed  boundar 
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Boundary 
Condition 
Along  this  line 

X=0/-^ 

L=X  T=0 


/ 

Figure  1:  Initial  Condition  Prescribe 

Satisfying  the  prescribed  boundary  conditions  in  general,  the  study  of  pressure  waves  in  a fluid  propagation  of 
heat  and  unsteady  state  problems  lead  to  parabolic  type  of  equations.  In  this  paper  we  are  studying  the  finite  difference 
algorithm  for  solving  parabolic  partial  differential  equation. 

1.2  Two  Dimensional  Heat  Equation 

The  heat-conduction  equation  can  be  applied  to  more  than  one  spatial  dimension.  For  two  dimensions,  its  form  is: 


Open  ended  domain  R 


Boundary  along 
(his  line 


d2u 

.dx2  dy2/ 


Where,  ff 


GZt 


Example:  1 (Solve  One  Dimensional  Heat  Equation) 

Find  the  analytical  solution  of  the  parabolic  equation 


= 2 ut  Whpn  u(0,  t)  = u( 4, 0 = 0 anH  u(x,  0)  = x(4  - x),  taking/f  = 1 and  c 


2 _ 


1 

2 


Find  the  Values  up  to 


t = 


5 


Solution 


It  is  already  given  parabolic  equation 


U 


= 2Ut 


1 d2u  Lili 

2 dx:2  lit- 


and  also  given  the  boundary  condition. 


lt(G,£)  = u(4,t)  — 0 


the  initial  condition 


i l(x,  0)  = :v(4  — x), 


we  can  write  parabolic  equation  in  the  form. 


dw  1 

tit  2 dx2 


C 


2 


Where 


1 

2 


0) 


Impact  Factor  (JCC):  4.2949 
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Let 


u=  XT 


X x T”  f 

Where  is  the  function  of  1 only  and  is  the  function  of  only 


Differentiating  (2)  partially  W‘  r‘  ^ ^ we  get 


J^rfT,r 

3t 


dzu 
dx 2 


= X"T 


„2 


Let  (3)  and  (4)  equal  P 

XT'  =c2  XUT  = -p2 

Separating  variables,  we  get 
2 


x”  _ i n 
X ~ c*  T 


* - = -P 


Xr'  = and  T = — 

dx2  dt 


Where 

r = -c2  p2r 


dT  j 2t 

— = — C VI 

dt  ^ 


DT  = - c 1 p2T 


m = —c1  p2 


T= 

Now, 

X”  = - p 2 X 

d7  X 2 v 

— - = -pz  X 

dxz 

D2X  = -p2  X 
m = ± ip 

X = c2cospx  + c3sinpx 


(2) 


(3) 

(4) 
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c = - 

We  have,  2,  put  the  value  of  A and  ‘ in  equation  (2)  we  get. 


V 7*  . 

“ onrl  ■* 


it  = (c2  cospx  + cs  sinpr  ) cte  z 


u(0,t)  — 0 


x = 0 

, where  in  equation  (5),  we  get 


Now,  put  the  value  of 
~Vzt 

0 = C 2 ^ 

c = o,  , t t , r=  c1e~c“p2t  c1  ^ 0, 

Let  2 ’ arbitrary  constant  we  have  i so  1 


it  = (c3  sinpx}  cte  z 


Put  the  value  of 


uC4,t)  = 0when.x  = 4, 


in  equation  (6),  we  get 


-ft 


0 =c3  sin 4 pc±e  z 
3in  4p  = 0 = sinriTt 

nw 

V = — 
r 4 

Put  the  value  of  p in  equation  (6) 

M=  c3sm{- }x  z 


(n  air% 

m = c3  stn(^)x  e’  32  f 


Where  C3  ~ 2 ^ 

/trBgBY 

u=  e V 3-  ■ ' 

#m 

Now  we  find  the  value  of  by  Fourier  series. 


K = 7 J0  /(*)  sin(“)  dx 

2 , jucr,  , 

= - L x(4  — x)  smf — dx 

4 ' * * 4 ' 

, 1 f4,  . /rmx^  , 1 r4  -j  , , 

K =7  J0,  4rsin(— )dx-  - r sm(— ] dx 


(5) 


(6) 


(7) 


(8) 


(9) 


Impact  Factor  (JCC):  4.2949 


Index  Copernicus  Value  (ICY):  3.0 


Numerical  Methods  for  Parabolic  Paritial  Differential  Equations 


103 


Integration  first  part  from  equation  (9)  by  the  form  of  integral, 

4 


h = 2 


— x- 


sin 

c 

Ban 

4 J 

2 

1 

L = -32- 


(10) 


Integration  second  part  from  equation  (9) 


/3  “ 2 


— X 


2 COs{—  ' I 

Of)  J 


— i-x 

mr,  I 


sin 

1+  S 

\ 4,  J 

( 

n tf\ 

4 J 

f 

]i  mr 

\ 4 / 


cos  me  . 8 I cosrwr  1 
-16  -755T  + - 


t'js  Hi  i o i N.v’art.'v 

IfT 


h = 


C0STCJT  , 64  , 1 

-32  + - (COSTIU  — 1 ) 

w r narra  ' J 


(11) 


Putting  the  value  of  (10)  and  (1 1)  in  (9)  we  get, 
i nnEitfcmt  , uuariw  64 

bn  = -32— 1-  32  —{costlk  — 1) 


b„  = — (cos  me  — 1) 


When  U = 0 


b0  = 0 


128 

= JTa 

Put  value  of  ^1  in  equation  (8),  we  get 


128  . -f— V 

u = — sm(—)  xc  1:1  * 


(12) 


Where 


n — 1 


4,123196  ^ 

It,  = T= — x e 


-©■ 


ut  = 


V2 

4,128196  X 0,731615 


V 2 


= 2.1356 
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u?  = 4.128196  X e («)* 
u2  = 3,032585 

u3  = 4.128196  e”Wf 

u2  = 2.1356 


Analytical  solution  for  ^ ' 

NUMERICAL  SOLUTION 

Now  We  Solve  (Example:  1)  by  Numerical  Methods 

Find  the  solution  of  the  parabolic  equation  2ut  wken  0 0 and 

1.1  (x,  0)  = :\'(4  — X),  taking  h = 1 Find  the  vajues  Up  to  r = 5 

Solution 


Given  a parabolic  equation  ^xx  ^ 

1 d2u  du 


2dx2  dt 

Also  given  the  boundary  condition,  ^ — $ — ^ 

and  the  initial  condition  u(xr  = x(4  — X), 

X = 1 2 34 

when  we  put  the  value  of  ' ' ’ \n  the  relation  of  the  initial  condition,  it  gives  that, 

u(G,0)  = 3,  <1,0)  = 4,  <2,Q)  = 3,  u(3,G}  = 3,  u(4,0)  = 0, 
h = 1,  k = 1,  and  c = * 

we  have  ‘ using  these  values  we  can  find 


a — 


fee2  1 


ft2  2 By  Schmidt  method  a must  be  z 


By  Schmidt  Formula  the  following  Results  are  Obtained, 
When  ! = 1j2,3  a|)d  } = 0, 


T-h  = - [0  + 4]  = 2 
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^=i[3  + 3]  = 3 

u2  = 2 = 2 
When  i = 1'2'3  andj  = 3- 

u5  = i[Q  + 3]  = 1.5 

z 

u6=  i[2  + 2]  = 2 

2 

«7=  \ P + 0]|  = 1.5 

Li 

When  1 = 1'2'3  aUdi  = 2 

lr 

tig  — — [0  + 2]  = 1 
2 

«!,=  ^ [1.5 + 1.5]  = 1.5 

Jiu-  j[2  + 0]-l 


Now  Solve  by  Crank-Nicolson  Formula 

— aut_  jj+j  + (2  + 2 a)uij+  ± — &ui+1/j+1 


Wi-ij  + (2  - 'ZaJiiij  + aui+liJ 


a = 


Put 


1 

2 


Ui-Xj+1  ^ ui,j+l  Mi+l,j+l  ui-l,j  2 uij  ui+l,j 
When  1 = 0^,3  ClTldj  = 0, 


— 0+6  ut/1  ~ U21  = U0  + 2 Uj  + U2 

- Ultl  + 6 u2/1  - u2rl  = ut  + 2 u2+u2 
-U2/1  + 6 U2/1  - u41  = u2  + 2 u2  + u4 


When  the  initial  conditions  are  given  in  ® ' 00'  O^X  we  get. 


CO 

00 

(iti) 
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6 itj  — u2  + Gi£3  = 10 

— ut  + 6u2  — u2=  14 

— u2  + 6 u3  + 0 u4  = 10 

When  we  solve  the  equations  (H)1  (Hi)>  by  manual  or  by  calculator  it  gives  that, 

u3  = 2.176470588,  u2  = 3.058823529,  u2  = 2.176470588 

When  ^ 1^4^  CUldj  1, 

In  the  same  way  can  we  find, 

u4  = 1.615916955,  u5  = 2.283737023,  u6  = 1.615916955 

RESULTS 


CO 

(“3 

(Hi) 


Table  1:  First  Step  Using  the  Analytical  Solution  and  Crank-Nicolson, 

Bender-Schmidt 


Analytical  result 

2.1356 

3.032585 

2.1356 

Schmidt 

2 

3 

2 

Crank 

2.176470588 

3.0588233529 

2.176470588 

Table  2:  Second  Step  Using  the  Crank-Nicolson  and  Bender-Schmidt 


Schmidt 

1.5 

2 

1.5 

Crank 

1.615916955 

2.283737023 

1.615916955 

Ut 

u2 

CONCLUSIONS 


On  the  basis  of  the  above  discussion  we  get  the  result  obtained  by  analytical  methods  is  always  providing  accurate 
solution  but  numerical  solution  always  providing  approximate  result.  But  among  these  numerical  methods  Crank-Nicolsan 
method  was  providing  fast  convergence  in  comparison  to  Bender-Schdimit  method.  Since  it  is  not  possible  to  solve  every 
partial  differential  equation  analytically  so  numerical  methods  providing  a good  agreement  in  those  cases  where  solutions 
not  exist  or  we  are  unable  to  solve  partial  differential  equation  analytically. 
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